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Abstract. In this Comment we compute the contributions of the radiation reaction 
force in the 2.5 post-Newtonian (PN) gravitational wave polarizations for compact 
binaries in circular orbits, (i) We point out and correct an inconsistency in the 
derivation of [1]. (ii) We prove that all contributions from radiation reaction in the 
2.5PN waveform are actually negligible since they can be absorbed into a modification 
of the orbital phase at the 5PN order. 
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1. Introduction 

The second and a half post-Newtonian (2.5PN~ c~^) waveform for inspiralling compact 
binaries moving in quasi-circular orbits was computed by Arun et al [1]. Starting from 
the expressions of the radiative multipole moments given by Eqs. (3.4)-(3.6) of [1], 
and of the source moments given by (3.16)-(3.18) of [Ij, the wave form is made of the 
instantaneous terms (Eqs. (5.1)-(5.4) of [Ij), the hereditary memory-type contributions 
(Eqs. (4.24) of and the tail contributions (Eqs. (4.38) corrected by the published 
Erratum [Ij). In this Comment we are concerned with the piece of the instantaneous 
waveform that is the contribution from radiation reaction (RR) at 2.5PN order in the 
dynamics of the inspiralling compact binary. We point out and correct an inconsistency 
in the derivation of RR terms in [1] , but argue that in fact the RR terms are negligible 
in the 2.5PN waveform, since they can be absorbed into a 5PN order contribution to 
the orbital phase evolution. 



Radiation reaction in the 2. 5PN waveform 



2 



2. Quasi-circular inspiral at 2.5PN order 

Adopting the conventions of [1], r = |x| is the binary separation with x = yi — y2 
the vectorial separation between the particles, v = vi — V2 is the relative velocity, 
m = mi + 7712 is the total mass of the binary system (distinct from the mass monopole 
moment M given by the ADM mass of the system), and u = mim2/m^ is the reduced 
mass divided by the total mass. 

In modeling the orbital motion of the binary at 2.5PN order, Ref. [1] stressed 
the importance of including the radiation-reaction force in the 2.5PN expression of the 
binary acceleration. Consider an arbitrary orbit confined to the x-y plane. The relative 
position, velocity, and acceleration are given by 

X = rn, (1) 
V = fn + ruX, (2) 
a = (r — rco'^)n + {rcu + 2ruj)X. (3) 

Here A = z x n where z is the unit vector along the z-direction. The orbital frequency 
u is related as usual to the orbital phase via u = (p. 

Through 2PN order, it is possible to model the motion of the binary as a circular 
orbit with the solution r = r = u = and ru^ = — n ■ a. Detailed calculations of the 
2PN equations of motion for a circular orbit (in harmonic coordinates) yield 

At 2.5PN order, however, the effect of the inspiral motion must be taken into 
account. The leading order contribution to the inspiral can be obtained by examining 
the Newtonian orbital energy, 

E = -i^, (5) 
2 r ' ^ ^ 

and the leading order gravitational-wave luminosity, 
dE _ 32 G^m^z/^ 

~ ~dt ~ '5 r^c^ ^ ^ 

Here we assume that the energy radiated by the gravitational waves is balanced by the 
change in the orbital energy. This yields 

dE/dt 64GWz/ 



dE/dr 5 r^c^ 
Similarly the orbital frequency changes by (using Gm = r^u'^ at leading order) 



(7) 



dE/dt 9QG7niy /Gm\ 



5/2 



dE /duo 5 \ rc^ 

while the orbital phase = / ujdt, found by integrating ([8]), is 

1 fGTTl^ '^''^ 



) , (8) 



32z/ V rc2 ' ■ 
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Substituting ([7]) and into and noting that r ~ 0{c~^^) is of the order of 

the square of RR effects, the expressions for the 2.5PN inspiral relative velocity and 
acceleration in harmonic coordinates are obtained, 

^ = rujX-——^n, (10) 
2 32 G^m^u , , 

where lo is given by (jlj). As we shall detail the 2.5PN RR terms in both the inspiral 
velocity ffTOl) and acceleration ffTTl) should be substituted into the gravitational waveform 
at 2.5PN order. 

3. Computation of polarization waveforms 

In the leading quadrupole approximation the gravitational waveform is given by 

2G 



Im + 0{c-')\, (12) 

where R is the distance to the observer and Vijki is the transverse-traceless (TT) 
projection operator Vijki = VikPji — \VijVkh with Vij = 6ij — NiNj and N = (Ni) the 
radial direction from the source to the observer. Here Iki is the second time- derivative 
of some appropriate source quadrupole moment defined from a general post-Newtonian 
multipole moment formalism. The remainder 0{c~^) indicates higher PN corrections 
coming notably from all the higher multipolar orders. 

Defining an orthonormal triad (N, P, Q) where P and Q are unit polarization 
vectors transverse to the direction of propagation, the polarization waveforms h+ and 
hx are computed from the waveform f|T2l) by 

h+ = l{P,P,-Q,Q,)hj;', (13) 

hx = UP^Qj + PjQ^)hJ;'. (14) 

If the orbital plane is chosen to be the x-y plane with the orbital phase measuring the 
direction of the unit vector n = x/r along the relative separation vector, then 

n = cos 0x + sin 0y. (15) 

Ref. [1] has chosen the polarization vector P to lie along the x-axis and the observer to 
be in the y-z plane with 

N = siniy -|- cosiz, (16) 

where i is the orbit's inclination angle (0 < i < vr). With these definitions P lies along 
the intersection of the orbital plane with the plane of the sky in the direction of the 
ascending node, and the orbital phase is the angle between the ascending node and 
the direction of body one (say). The rotating orthonormal triad (n. A, z) describing the 
motion of the binary is then related to the polarization triad (N, P, Q) by 

n = cos(/)P + sin0(QQ + SiN), (17) 

A = - sin0P + cos0(qQ + SjN), (18) 

z = -SiQ + QN, (19) 
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where we pose Cj = cosi and Sj = sini. 



4. Radiation reaction contributions to the waveform 

All contributions arising from the gravitational RR are contained in the leading order 
quadrupolar waveform given by ( fT2l) . They have three different origins. The first type 
of RR term is issued directly from the expression of the (post-Newtonian) quadrupole 
moment lij at order 2.5PN, and reads (see Eq. (3.16a) of [1]) 

lii = vmx^'x^^ H h — ^ — x<V>. (20) 

7 ra' 

The first term is the usual expression of the Newtonian quadrupole moment (the brackets 
surrounding indices refer to the symmetric-trace- free projection). The dots indicate the 
IPN and 2PN conservative terms, and terms strictly higher than 2.5PN. Taking two 
time-derivatives of fl20l) we get 

1 no r^3™4, .2 

'L = 2um(v<V> + x<'a^>) + ^ , , x<V>. (21) 

The second type of RR term comes from inserting the expression of the acceleration 
a given by (fTTj) . which produces another RR term which modifies the term in (12T!) as 

Gm ,.\ 1408^3^^2 
a;<V> + — 

The sum of these RR contributions is exactly what has been computed in pQ, where it 
was denoted by p^j^^'* and given by Eq. (5.2). We now replace fl22l) into the waveform (fT2l) 
and compute the two GW polarizations according to (fT3l) - (fT^ . hence 

^i^j QiQj 

\ r / . . t .\ I ZJ.I 1^ / ^"m "7/- .1 

(23) 



4- = 2z/m(w<V> ^x<V> + — ^^x<V^>. 22 



h 2^ 



_ 2 



2 



o I i i Gm i A 1408GWz/2 ^ , 
/z/m( f f —X x-^ ] + ■ ■ ■ — — X 



For ease of notation, in Eq. ( 123|) above and similar equations later the first row (line) 
corresponds to + and the second row (line) to x. The RR term in (l22l) has already 
been included in the final result of [I] . Its contribution to the polarization waveforms is 
given by 

2G.^Jf(^ + ^^^)^^''^^^^^^' 
Suh^,^ = -^{ (24) 

' [ -^(2c.)z.x^/2cos20, 

where we pose x = {Gmuj / c'Y^^ . With the notation 5i2 we remind that this term was 
made of two distincts pieces. 

However, let us show there is also another contribution to RR that has been 
overlooked in [T]. Indeed the instantaneous terms in the waveform hj^ (Eqs. (5.1)- 
(5.4) of [1]) are given in terms of the relative position n and velocity v of the binary, 
as well as the PN parameter 7 = Gm/{rc^). In computing the polarization waveforms, 
i.e. projecting out hj^ to get /i+,x, Ref P has substituted, at the last stage of the 



Radiation reaction in the 2. 5PN waveform 



5 



computation, v = ruX in the waveform to obtain Eqs. (5.9)-(5.10) of [T]. This is correct 
for all the terms except the two leading order "Newtonian" terms in fl23l) for which one 
must use the true expression of the velocity ( |T0|) including RR inspiral. We find that 
substituting the inspiral velocity (fTOj) into the leading terms of the waveform (l23l) yields 
the following additional contribution to the polarization waveforms, 

64 



2Gvmx 



'l + cl)vx^/^ sin 20, 



64 

-y(2ci)z/x^/' cos; 



(25) 



Notice that both results ( |24l) and ( !25l) have the same structure. Since the 
contribution ( l25l) has not been taken into account in [1] we find that the following 
terms in the 2.5PN polarization waveforms of [1] (Eqs. (5)-(6) of the Erratum) are to 
be changed from 



H 



(2.5), 



I old 



+ sin 2il} 



+ Ci cos 2il) 



9 14 



7 



/96 8 



28 



+ 



22 . 
2 - -c? + . 



154 94 . 



(26) 



to 



H 



(2.5) 



+ sin 2il) 



■ + Ci cos 2iIj 



9 14 2 7 4 / 56 2 



+ 



28 



+ 



22 2 
2 - — c,^ + 



282 94 . 



(27) 



+ ■ 



5 ' V 5 

We recall that the phase variable ip differs from and is given by Eq. (5.6) of [1]. The 
difference between ip and is at order 4PN at least (see pTj for discussion). Note that the 
new correction (125|) is in terms that vanish in the limit u 0, and so is still consistent 
with the results of black hole perturbation theory [2j. 



5. Phase modulation due to radiation reaction 



We are now going to show that all the RR contributions are in fact negligible in the 
2.5PN waveform. From the sum of the previous results (p^ and ( [25|) we end up with 
the following total contribution due to RR in the GW polarizations: 

^(1 + q^)z/x^/^ sin; 
^^\2ci)ux''/^ cos 20. 
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Comparing this result with the GW polarizations at Newtonian order, 



h 



(N) 
+,x 



2Gh'mx 
c^R 



-(1 + c.)cos: 

-(2ci) sin 20, 



(28) 



(29) 
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we see that the RR terms can be absorbed into a redefinition of the phase variable as 



80 



5/2 



(30) 



This means that the 2.5PN waveform will take exactly the same expression as if we 
had neglected the RR contributions ( l28l) . but with phase variable $ in place of (p. Now 
the point is that the added term in fl30|) represents an extremely small modification of 
the phasing. Indeed, we take into account the expression for the phase as a function of 
frequency at the leading order {i.e. due to leading order RR effects), which has been 



computed in ([9]) and whose order ~ c' is the inverse of the order 
Thus we find that the redefined phase is equivalent to 

p 2560 



32z/ 



1 + - 
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of RR effects. 



(31) 



This means that the RR terms seen as modulations of the phase evolution, contribute to 
the phase much beyond the 2.5PN order, namely at order 5PN~ beyond the leading 
phase evolution. Indeed there is in principle no point in including these terms because 
they are comparable to 5PN terms in the phase that are unknown — only the phase 
evolution up to 3.5PN order is known. 

We conclude therefore that the RR terms can in fact be neglected in the 2.5PN 
waveform for circular orbits. This is similar to what happens in the energy flux 
for circular orbits where we know that the 2.5PN radiation reaction gives finally a 
contribution only at 5PN order [3] (the 2.5PN terms in the energy fiux are only due to 
tails). Hence we can by the redefinition of the phase variable (|30ll completely remove 
the RR contribution given by fl28|) . However we can also decide to include such terms 
in the 2.5PN waveform (as usual there are many different ways of presenting PN results 
at a given order of approximation). Here we simply propose to include the RR terms 
as they are in the templates of binary inspiral, i.e. by correcting the inconsistency 
in [Ij (which has by the above argument no physical effect at 2.5PN) and keeping the 
corrected waveform in the form given by Eqs. fl27j) . For completeness, we note that if 
we choose to include all the RR terms into the phase redefinition $, Eq. ( 1271) will then 
be modified into 



■ + sin 2^ 



H 



/(2.5) 



9 14 



7 



/64 408 



28 



+ ■ 



■ + Ci cos 2* 



\. 22 2 
2 - —c^ + y 



374 94 



(32) 



■3^ + Y^^ 



+ ■ 



where ^ is equal to ip plus the added RR contribution given by the second term in fl30|) . 
and where if) itself differs from by Eq. (5.6) in pLj. (All the other terms in the waveform 
are then to be expressed with the same phase variable ^ .) 
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